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ABSTRACT 

Context. The theory of stellar evolution can be more closely tested if we have the opportunity to measure new quantities. Nowadays, 
observations of galactic RR Lyr stars are available on a time baseline exceeding 100 years. Therefore, we can exploit the possibility 
of investigating period changes, continuing the pioneering work started by V. P. Tsesevich in 1969. 

Aims. We collected the available times of maximum brightness of the galactic RR Lyr stars in the GEOS RR Lyr database. Moreover, 
we also started new observational projects, including surveys with automated telescopes, to characterise the O-C diagrams better. 
Methods. The database we built has proved to be a very powerful tool for tracing the period variations through the ages. We analyzed 
123 stars showing a clear O-C pattern (constant, parabolic or erratic) by means of different least-squares methods. 
Results. Clear evidence of period increases or decreases at constant rates has been found, suggesting evolutionary effects. The median 
values arey8=+0.14 d Myr~' for the 27 stars showing a period increase and/?=-0.20 dMyr~' for the 21 stars showing a period decrease. 
The large number of RR Lyr stars showing a period decrease (i.e., blueward evolution) is a new and intriguing result. There is an excess 
of RR Lyr stars showing large, positive ft values. Moreover, the observed fi values are slightly larger than those predicted by theoretical 
models. 
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1. Introduction 

RR Lyr variables are low-mass stars in a core helium burning 
phase; they fill the part of the HR diagram where the horizon- 
tal branch intersects the classical instability strip. The crossing 
of the instability strip can take place in both directions; as a 
consequence, the periods will be either increasing, if the stars 
evolve from blue to red, or decreasing, if they evolve from red to 
blue. Despite its importance as a test for the stellar evolution the- 
ory, the observed rate of the period changes is still an unknown 
quantity. It has been measured in globular clusters (e.g., Smith & 
Sandage |l981| Lee l 19911 Jurcsik et al. 120011 1. comparing periods 
determined in different epochs. However, this task has not been 
undertaken yet on the wide population of galactic RR Lyr stars, 
though most of them have been studied for tens of years, several 
of them since the end of the XIX th century, and therefore the 
available data are almost continuous. One of the critical points 
is to determine the importance of the negative rates, i.e., the pe- 
riod decreases. Such rates have been observed in some cluster 
stars (Jurcsik et al. 2 001 \, but it is unclear in which evolutionary 
phase they occur (Sweigart & Renzini [l9791 >. 



In the present paper, we derive the period variation rate of 
the best observed RR Lyr stars belonging to the RRab sub- 
class. For this purpose we use the GEOS RR Lyr database 
which is described in Sect. [2] The observations of the TAROT 
telescopes coordinated in the GEOS RR Lyr Survey give 
a strong impulse to this analysis, providing a large num- 
ber of times of maximum brightness in the recent years. 
They complete the effort of the amateur observers, in par- 
ticular in the European associations BAV (Bundesdeutsche 
Arbeitsgemeinschaft fiir Veranderliche Sterne) and GEOS 
(Groupe Europeen d' Observations Stellaires), which have been 
surveying these stars for years. These new observations are de- 
scribed in Sect. 12.21 Next, Sect. [3] shows how the data were anal- 
ysed. In Sect. |4j we analyse the stars with constant period vari- 
ation. Section [5] is devoted to the study of some particularities 
encountered during the present work: Blazhko and light-time 
effects. Finally, a discussion of the results is developed in Sect. [6] 
in the light of stellar evolution. 
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2. The GEOS RR Lyr database and survey: a V.P. 
Tsesevich's heritage 

In his book "RR Lyrae stars", V. P. Tsesevich d 19691 ) made a 
summary of the period variation behavior of RR Lyr stars from 
the observations made during the 60 years or so after the dis- 
covery of these stars. Nowadays, 40 more years have elapsed, 
allowing 100 years of data for some stars. It is timely to investi- 
gate what the data accumulated since Tsesevich's work bring to 
the understanding of the period variation of RR Lyr stars and of 
the stellar evolution of horizontal branch stars. 



2.1. The GEOS RR Lyr Database 

To complete Tsesevich's data collection, the ama- 
teur/professional association GEOS has built a database 
aiming to put together all possible RR Lyr light maximum times 
published in the literature. The publications from the end of 
the XIX th century up to today were scanned for this purpose. 
The database continues to be fed by recent observations from 
amateur astronomers of the European groups GEOS and BAV. 
To date, the databas43 contains about 50000 maximum times 
from more than 3000 RR Lyr stars. 

The stars concerned are galactic RRab and RRc. The col- 
lected times of light maxima were obtained from observations 
using various techniques, either visually, with electronic devices, 
or photographically, and by hundreds of observers. This hetero- 
geneity is taken into account in the analysis of the O-C diagrams 
which are the basis of the present paper. The O-C diagrams are 
obtained plotting the differences between the observed and cal- 
culated times of maximum brightness (i.e., the O-C values) vs. 
the elapsed cycles. Multicolor photometry of RR Lyr stars (e.g., 
the extensive UBV light curves obtained by Olah & Szeidl [1978l 
and Szeidl et al. 1 19861 shows that the colour effects on the deter- 
mination of the times of maximum brightness are negligible with 
respect to the period variations we are searching for, since they 
are much smaller than the errors on the determinations them- 
selves. It is also obvious that the precision of modern observa- 
tions, either with photoelectric detectors or CCDs, is much better 
than visual and photographic observations. However, the uncer- 
tainties of the visual and photographic determinations are much 
less than the O-C variations we are studying. A particularity of 
photographic observations is that the corresponding published 
maximum is often not the result of the analysis of a lightcurve 
obtained during one night, but rather one single "bright" point. 
These observations were useful at the time when people searched 
for the pulsation period of the star, but unavoidably this method 
introduces an additional scatter in the O-C diagrams. 

2.2. The TAROT Telescopes 

The telescopes TAROT (Telescope a Action Rapide pour les 
Objets Transitoires — Rapid Action Telescope for Transient 
Objects; Bringer et al. 1 19991 1 are automatic, autonomous obser- 
vatories composed of 25-cm telescopes whose first objective is 
the real-time detection of the optical transient counterparts of 



1 The GEOS database is freely accessible on the internet at the ad- 
dress http://dbrr.ast.obs-mip.fr/ This site is hosted by the Laboratoire 
d'Astrophysique de Toulouse-Tarbes, Observatoire Midi-Pyrenees, 
Toulouse, France 



GRBs in the visible (Klotz et al. 12006 PI. One of these telescopes 
is located in the northern hemisphere at Calern Observatory 
(Observatoire de la Cote d'Azur, University of Nice, France), 
and the second one is in the southern hemisphere at ESO, La 
Silla, Chile. The telescopes have a field of view of 1.86x1.86 
deg 2 . During operation no human action is necessary, neither for 
calibration and science frames nor for data processing. The roof 
of the buildings are retractable to allow rapid movement to any 
direction in the sky. TAROT telescopes are able to point in any 
direction in less than 5 sec (including CCD operations and fil- 
ter movements). They can detect stars of magnitude V=16.2 in 
lOsec exposure time, or V=18.2 in 60 sec. TAROTs are able to 
monitor a GRB 8 sec after being triggered, about 20 sec after the 
gamma ray detection (Boer et al. l2.0011 >. 

A dedicated software program optimises the process: it is 
able to decide at any time what is the best observation to do next, 
given program objects, priorities, alerts and constraints from 
the astronomer such as measurement frequency (Bringer et al., 
|1999l >. Weather stations are linked to the TAROT computers. If it 
rains or if the temperature reaches the dew temperature or the hu- 
midity becomes higher than 90%, the roof remains closed. The 
system performs an automatic data reduction (dark and flat cor- 
rections, astrometric and photometric calibrations) in real time 
(see Damerdji et al. 2007). In the case of the images dedicated 
to the RR Lyr survey, a catalogue of the measurements is up- 
dated after each observation and is available on the web about 
2 min after the image has been taken. 

2.3. The GEOS RR Lyr survey with TAROT telescopes 

In addition to the detection and localisation of GRBs, TAROTs 
follow the CNES artificial satellites, and a fraction of the observ- 
ing time is left to other astronomical programs. Therefore, the ef- 
forts made by hundreds of observers over one century to monitor 
the maximum times of RR Lyr stars may be continued nowadays 
by much more efficient techniques, thanks to the use of these 
small robotic telescopes. Indeed, the authors started a monitoring 
of maxima of RRab stars, called GEOS RR Lyr survey, on 2004 
January 29 using the northern TAROT and in September 2006 
using the southern one. The capabilities of the system and the 
time-sharing allow us to schedule the observation of up to 6 max- 
ima of RR Lyr stars per night and per telescope. The schedule is 
sent to the control sofware once per month. The aim is to obtain 
at least one maximum per star with V < 12.5 (Calern observa- 
tory) or with V < 13.5 (Chile) and per year, in order to check 
the period stability on time scales of decades. In practice, the 
number of maxima observed for each star of the programme is 
typically 5-10 per year. Only stars from the General Catalogue 
of Variable Stars (GCVS, Kholopov et al. 11988b have been ob- 
served so far: this catalogue contains stars which have often been 
observed for more than 50 years. In 3 years of operation, from 
January 2004 to December 2006, TAROT-Calern observed 1288 
maxima on 130 stars; TAROT-Chile observed 556 maxima on 
126 stars during the first 6 months of operation. 

For each maximum, the telescopes obtain an average of 50 
measurements in 4.8 hours. The number of frames having con- 
tributed to determine the maximum times is 98000 in three years. 
The times of maximum brightness are determined by fitting the 
data points by means of a cubic spline function with a non-zero 
smoothing parameter. The smoothing parameter was selected 



2 TAROT telescopes were funded by the Centre National de la 
Recherche Scientifique (CNRS). The project is led by a team of French 
laboratories (CESR, LATT and OHP). 
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case-by-case, depending on the number of points and their scat- 
ter. The error bars are estimated from the data sampling. The 
nominal sampling is composed of two consecutive 30-sec expo- 
sures taken every 10 min in a time interval of 2 hours, centred 
around the predicted maximum time and 20 min outside this in- 
terval. The resulting mean separation between 2 measurements is 
5.8 min. However, the nominal sampling may be altered by local 
events (e.g., weather or telescope operation). For a well covered 
maximum, the mean error bar is about 0.003 d (4.3 minutes); 
usually, error bars range from 0.002 to 0.010 d. 

2.4. Survey statistics 

In order to check the efficiency of the survey, we analysed the 
production of maximum times during the first three years of 
TAROT-Calern. The RR Lyr stars in the GCVS with V > 9.0 
(at minimum) and with 6 > -10° are sorted in bins of 0.5 mag 
width. RR Lyr itself is not considered, since this bright variable 
(7.06-8. 12 V) is currently monitored with several techniques (vi- 
sual, photographic, photoelectric and CCD photometry, medium 
and high resolution spectroscopy). Since the magnitudes at mini- 
mum come from the GCVS, the photometric systems are hetero- 
geneous, but close to Johnson V magnitudes. We also compare 
the statistics of the TAROT (GEOS) production (from 2004 to 
2006) with those of the observers of BAV group (from 2002 
to 2005). Both groups published their maximum times in the 
Information Bulletin on Variable Stars: see Le Borgne et al. 
I2D041 l2005al I2005bl l2006al I2006bl and 120071 for GEOS max- 
ima, and Ager er and Hiibs cher, [20021 [20031 Hiibscher [20051 
Hiibscher et al. l20031l2006l for BAV maxima. BAV obtains this 
result using several telescopes operated by about 20 observers. 

Tableland Fig. Q] illustrate the results of the statistical anal- 
ysis. Indeed, the results show that the two groups put together 
have a completeness of 97% up to V=12.5 (magnitude at min- 
imum); 93 of the 97 RRab stars listed in the GCVS have been 
observed. Individually, the GEOS and BAV surveys have a com- 
pleteness of 88% and 81%, respectively. We note that, for the 
GEOS survey with TAROT telescopes, V=12.5 is a natural limit 
for a 25-cm telescope. Moreover, both surveys have a decrease 
in completeness beyond this magnitude imposed by the typical 
diameter of the telescopes (20-30 cm mirror diameter), while 
a second maximum occurs in the BAV histogram at magnitude 
13.5, corresponding to the use of larger telescopes by some ob- 
servers. 

3. Data analysis 

We extracted from the GEOS database the RRab stars for which 
the times of maximum brightness span more than 50 years. 
Actually, the observations of RR Lyr stars started at the end 
of the XIX th century and for several objects the time base- 
line exceeds 100 years, with data in three different centuries. 
Throughout this time the observing techniques were very differ- 
ent (visual, photographic, photoelectric, CCD) and for each tech- 
nique the accuracy of the data was not homogeneous, depend- 
ing on the observer, the equipment and the number of measure- 
ments. Unfortunately, all that information was seldom known 
in the case of very old observations (often disseminated in sev- 
eral papers), which in turn are very important for establishing 
the O-C trends. In these conditions, the attribution of weights to 
the maxima observed with different techniques was not practi- 
cable in a rigorous way. Moreover, when viable, this weighted 
procedure yielded results not significantly different from those 
obtained without weights; indeed, in these cases the datasets 



Table 1. Number of RR Lyr stars per magnitude (at minimum) 
bin observed from 2004 to 2006. 



Magnitude Number of stars 

at minimum 

in GCVS BAV+GEOS GEOS BAV 

2002-2006 2004-2006 2002-2005 



9.5 - 10.0 


1 


1 


100% 


1 


1 


10.0 - 10.5 


9 


8 


88% 


8 


7 


10.5 - 11.0 


9 


9 


100% 


9 


8 


11.0- 11.5 


17 


17 


100% 


13 


15 


11.5 - 12.0 


23 


21 


91% 


19 


17 


12.0- 12.5 


38 


37 


97% 


35 


31 


12.5 - 13.0 


38 


29 


76% 


21 


22 


13.0- 13.5 


59 


46 


77% 


19 


41 


13.5 - 14.0 


69 


34 


49% 


4 


33 


14.0 - 14.5 


93 


23 


24% 


3 


21 


14.5 - 15.0 


108 


24 


22% 


2 


22 


15.0- 15.5 


164 


20 


12% 


2 


18 


15.5 - 16.0 


188 


13 


6% 


3 


10 


16.0- 16.5 


222 


5 


2% 





5 


16.5 - 17.0 


187 


2 


1% 





2 




magnitude at minimum 

Fig. 1. Number of RR Lyr stars per magnitude (at minimum) bin ob- 
served from 2004 to 2006. The observing material shown by the light 
grey histogram in the main window of the figure includes the GEOS 
RR Lyr Survey and the BAV CCD observations published in the IB VS. 
The dark grey histogram represents all the RR Lyr stars included in the 
GCVS. The enclosed histogram shows details of the contribution of the 
GEOS RR Lyr survey (filled blue) and the BAV (hatched blue). 



Table 2. Inventory of the O-C patterns in the RRab stars of our 
sample 



O-C pattern 


N 


Representative cases 


Constant period 


54 




Increasing period 


27 


BN Aqr, RR Leo 


Decreasing period 


21 


SW And, AH Cam, V964 Ori 


Irregular slow variations 


17 


UY Boo, RW Dra 


Single abrupt change 


4 


SZHya 


Total 


123 
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Table 3. Refined linear elements for stars showing a constant period. The note B beside the value of the standard deviation indicates 
the stars having a known Blazhko effect. 



Star N max Time coverage Epoch Period s.d. 

[years] [HJD - 2400000] [d] [d] 



AT 


And 


66 


101 


(1906 


2007) 


36109.5000 


± 0.0030 


f\ f 1 /"A 1 A A C 

0.61691445 


± 0.00000014 


0.0169 


DE 


And 


29 


64 


(1935 


■1999) 


39764.5234 


± 0.0033 


0.45363536 


± 0.00000019 


0.0176 


GV 


And 


54 


62 


(1937 


■1999) 


40171.3320 


± 0.0050 


0.52809310 


± 0.00000033 


0.0319 


OV 


And 


100 


75 


(1929 


2004) 


39469.2891 


± 0.0012 


f\ A ^ f\ C {A 1 1 A 

0.47058114 


+ 0.00000007 


0.0111 


TZ 


Aqr 


34 


92 


(1915 


2007) 


T-rr .11") / \ / v~> "1 

37548.9023 


± 0.0022 


0.571 19429 


+ 0.00000010 


0.0128 


YZ 


Aqr 


5 1 


1 12 


(1894 


2006) 


33343.2617 


r\ aai "> 

± 0.0033 


0.55193269 


± 0.00000017 


0.0218 


DN 


Aqr 


35 


71 


(1936 


2007) 


40915.3125 


± 0.0050 


0.63375407 


/ \ AAAAAATA 

+ 0.00000039 


0.0268 


TZ 


Aur 


140 


94 


(1913 


2007) 


37946.4922 


± 0.0006 


0.39167482 


± 0.00000002 


0.005 1 


BH 


Aur 


41 


100 


(1906 


2006) 


37352.4727 


± 0.0024 


0.45608902 


+ 0.00000007 


0.0086 


SS 


Cnc 


104 


96 


(1909 


2005) 


36201.3906 


± 0.0012 


0.36733848 


± 0.00000003 


0.0117 


RX 


CVn 


56 


98 


(1907 


2005) 


36904.2266 


± 0.0024 


0.54002589 


± 0.00000012 


0.0176 


AL 


CMi 


41 


82 


(1925 


2007) 


36607.3047 


± 0.0046 


0.55051273 


± 0.00000020 


0.0243 


UU 


Cet 


30 


82 


(1917 


1999) 


37582.4336 


± 0.0077 


0.60607427 


± 0.00000056 


0.0380 B 


SLJ 


Col 


25 


65 


(1937 


2002) 


39410.5234 


± 0.0064 


0.48735806 


± 0.00000035 


0.0302 


ST 


Com 


62 


90 


(1916 


2006) 


37652.6680 


± 0.0022 


0.59892809 


± 0.0000001 1 


0.0163 


W 


Crl 


23 


70 


(1935 


2005) 


39910.8633 


r\ r\f\ 1 a 

± 0.0014 


0.41201425 


± 0.00000007 


0.0063 


UY 


Cyg 


106 


105 


(1901 


2006) 


34862.3203 


± 0.0014 


0.56070578 


± 0.00000006 


0.0116 


BV 


Del 


28 


77 


(1928 


2005) 


39354.2539 


± 0.0037 


0.42345089 


+ 0.00000015 


0.0184 


BT 


Dra 


63 


101 


(1905 


2006) 


35868.4648 


± 0.0028 


0.58867335 


± 0.00000014 


0.0205 


SZ 


Gem 


74 


85 


(1922 


2007) 


38360.8438 


± 0.00 15 


0.50113541 


+ 0.00000007 


0.0108 


GI 


Gem 


63 


66 


(1941 


2007) 


36971.7148 


± 0.0016 


0.43326673 ± 0.00000007 


0.0122 


TW 


Her 


226 


94 


(1912 


2006) 


36784.3828 


/ \ f\f\f\ A 

± 0.0004 


0.39960003 


± 0.00000002 


0.0062 


BD 


Her 


62 


97 


(1900 


1997) 


33064.4766 


f\ AA^I 

± 0.0043 


0.47390789 


+ 0.00000025 


0.0317 


EE 


Her 


33 


57 


(1936 


1993) 


38853.6250 


± 0.0026 


0.49553666 


± 0.00000025 


0.0129 


V524 


Her 


122 


56 


(1930 


1986) 


36318.5156 


± 0.0032 


0.48186663 


± 0.00000022 


0.0326 


GO 


Hya 


39 


78 


(1929 


2007) 


35821.6641 


± 0.0052 


0.63643533 


± 0.00000028 


0.0317 


RX 


Leo 


23 


72 


(1935 


2007) 


38431.9336 


± 0.0055 


0.65341568 


± 0.00000027 


0.0115 


ST 


Leo 


96 


80 


(1926 


2006) 


39260.4219 


± 0.0007 


0.47798404 


+ 0.00000004 


0.0068 


WW 


Leo 


40 


72 


(1935 


2007) 


39801.6211 


± 0.0031 


0.60284597 ± 0.00000021 


0.0197 


AA 


Leo 


22 


52 


(1953 


2005) 


47271.3828 


± 0.0021 


0.59865355 


+ 0.00000023 


0.0088 


AX 


Leo 


35 


79 


(1927 


2006) 


37026.5742 


f\ f\f\ A 1 

± 0.0041 


0.72682691 


± 0.00000023 


0.0230 


V 


LMi 


34 


95 


(1912 


2007) 


36630.6094 


± 0.0018 


0.54391927 


+ 0.00000006 


0.0065 


VY 


Lib 


21 


77 


(1914 


1991) 


33773.3203 


± 0.0035 


0.53394097 


± 0.00000028 


0.0162 


TT 


Lyn 


50 


64 


(1943 


2007) 


43580.4023 


± 0.0032 


0.59743237 


± 0.00000023 


0.0224 


TW 


Lyn 


70 


48 


(1956 


2004) 


45022.4961 


± 0.0022 


0.48186052 


±0.00000013 


0.0187 


CN 


Lyr 


72 


105 


(1901 


2006) 


36079.3242 


± 0.0018 


0.41138276 


± 0.00000005 


0.0104 


IO 


Lyr 


95 


97 


(1909 


2006) 


36276.2305 


± 0.0013 


0.57712215 


± 0.00000007 


0.0113 


KX 


Lyr 


62 


65 


(1940 


2005) 


42019.2461 


± 0.0036 


0.44090429 


± 0.00000017 


0.0188 


LX 


Lyr 


33 


64 


(1940 


2004) 


36287.2930 


± 0.0026 


0.54548401 


+ 0.00000022 


0.0148 


RV 


Oct 


140 


91 


(1900 


1991) 


31874.5703 


± 0.0029 


0.57116264 


± 0.00000022 


0.0347 


ST 


Oph 


29 


98 


(1908 


2006) 


36079.3438 


± 0.0028 


0.45035613 


+ 0.00000011 


0.0144 


V531 


Oph 


25 


64 


(1928 


1992) 


39618.5156 


± 0.0033 


0.55365324 


± 0.00000022 


0.0166 


V2033 


Oph 


24 


52 


(1940 


1992) 


39270.5352 


± 0.0047 


0.56582946 


+ 0.00000043 


0.0201 


AO 


Peg 


23 


80 


(1925 


2005) 


37909.3984 


± 0.0021 


0.54724342 


+ 0.00000010 


0.0092 


BF 


Peg 


24 


73 


(1931 


2004) 


39369.1758 


± 0.0062 


0.49580538 


± 0.00000029 


0.0301 


CG 


Peg 


109 


51 


(1956 


2007) 


45576.4414 


± 0.0015 


0.46713728 


± 0.00000010 


0.0144 


DZ 


Peg 


63 


97 


(1910 


2007) 


36511.4062 


± 0.0043 


0.60734892 


± 0.00000026 


0.0313 


VY 


Ser 


41 


74 


(1932 


2006) 


39615.9297 


± 0.0049 


0.71409613 +0.00000038 


0.0225 


AN 


Ser 


106 


107 


(1899 


2006) 


34265.2461 


± 0.0014 


0.52207130 


+ 0.00000007 


0.0131" 


DF 


Ser 


42 


105 


(1900 


2005) 


34195.6641 


± 0.0019 


0.43779647 


+ 0.00000008 


0.0125 


UU 


Vir 


23 


97 


(1903 


2000) 


34509.2617 


± 0.0012 


0.47560632 


± 0.00000005 


0.0059 


WY 


Vir 


37 


87 


(1914 


2001) 


38060.6328 


± 0.0068 


0.60935342 


± 0.00000058 


0.0409 B 


AS 


Vir 


19 


84 


(1916 


2000) 


37052.4844 


± 0.0123 


0.55342495 


± 0.00000075 


0.0407 


AV 


Vir 


25 


93 


(1913 


2006) 


37016.4492 


± 0.0025 


0.65690929 


± 0.00000012 


0.0097 



were mostly composed of data from a small number of observers 
(often using the same technique) and the resulting differences 
among weights were small. In a later stage, we also tried to as- 
sign weights to different subsets considering their scatter from a 
fitting O-C curve (see Sect. |4|, and this procedure also yielded 
similar results. Therefore, we are confident that our conclusions 



are not affected by the constraint of having assigned the same 
weight to all the maxima. 

Particular care has been taken in the correct evaluation of 
the cycle values to be assigned to old maxima or after large 
gaps. Different O-C plots were obtained from different cycle 
countings, and abrupt jumps have been considered as suspicious. 
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Table 4. Parabolic elements for RRab stars showing a well-defined linearly increasing period. The note B beside the value of the 
standard deviation indicates the stars having a known Blazhko effect and discussed in the text. 



Star 




Time coverage 


Epoch 


Period 


Quad, term 


s.d. 


dP/dt 


P 


a 






[years] 


[rU U-Z^fUUUUUJ 


Ml 


[tu -aj 


Lei J 


rin-10 H/rll 
[1U ci/aj 


[a iviyi j 


|Myr J 


XX And 


182 


105 


54106.3508 


0.72276008 


1.427 


0.0167 


3.95 


0.144 


0.200 






(1902-2007) 


± 0.0027 


± 0.00000029 


± 0.066 




±0.18 


±0.007 


±0.009 


BN Aqr 


78 


110 


53979.4964 


0.46967415 


3.173 


0.0149 


13.51 


0.494 


1.051 






(1896-2006) 


± 0.0049 


± 0.00000026 


± 0.030 




± 0.13 


±0.005 


±0.010 


V341 Aql 


117 


94 


53999.4346 


0.57802322 


0.613 


0.0130 


2.12 


0.078 


0.134 






(1912-2006) 


± 0.0025 


± 0.00000027 


+ 0.056 




±0.19 


±0.007 


±0.012 


X Ari 


112 


93 


54107.2779 


0.65116811 


4.763 


0.0135 


14.63 


0.535 


0.821 






(1914-2007) 


± 0.0021 


± 0.00000028 


± 0.060 




± 0.18 


±0.007 


±0.010 


RS Boo 


325 


107 


54113.6372 


0.37734035 


0.268 


0.0093 B 


1.42 


0.052 


0.137 






(1900-2007) 


± 0.0011 


± 0.00000006 


± 0.007 




±0.04 


±0.001 


±0.004 


SWBoo 


52 


98 


53919.4937 


0.51355094 


4.639 


0.0080 B 


18.07 


0.660 


1.286 






(1908-2006) 


± 0.0019 


± 0.00000017 


± 0.030 




± 0.12 


±0.004 


±0.008 


UUBoo 


88 


99 


53904.4944 


0.45693375 


2.231 


0.0116 


9.77 


0.357 


0.781 






(1907-2006) 


± 0.0023 


± 0.00000017 


± 0.025 




± 0.11 


±0.004 


±0.009 


RWCnc 


95 


90 


53746.5042 


0.54721602 


1.937 


0.0169 


7.08 


0.259 


0.473 






(1916-2006) 


± 0.0035 


± 0.00000033 


± 0.054 




±0.20 


±0.007 


±0.013 


TTCnc 


97 


91 


54112.4014 


0.56345656 


1.229 


0.0132 s 


4.36 


0.159 


0.283 






(1916-2007) 


± 0.0023 


± 0.00000029 


+ 0.058 




± 0.21 


±0.008 


±0.013 


ANCnc 


29 


96 


53752.5259 


0.54316354 


0.555 


0.0122 


2.04 


0.075 


0.137 






(1910-2006) 


± 0.0043 


± 0.00000037 


± 0.067 




± 0.25 


±0.009 


±0.017 


SW CVn 


19 


95 


53478.5374 


0.44167101 


2.033 


0.0066 


9.21 


0.336 


0.762 






(1910-2005) 


± 0.0023 


± 0.00000024 


± 0.036 




± 0.16 


±0.006 


±0.013 


EZCep 


89 


105 


53750.3391 


0.37900204 


0.296 


0.0236 


1.56 


0.057 


0.151 






(1901-2006) 


± 0.0067 


± 0.00000027 


± 0.028 




±0.15 


±0.005 


±0.014 


RRCet 


113 


100 


54090.2963 


0.55302909 


0.231 


0.0075 


0.84 


0.031 


0.055 






(1906-2006) 


± 0.0014 


+ 0.00000011 


± 0.020 




±0.07 


± 0.003 


± 0.005 


TV CrB 


124 


104 


53903.4827 


0.58461776 


0.893 


0.0311 


3.06 


0.112 


0.191 






(1902-2006) 


± 0.0069 


± 0.00000047 


+ 0.074 




±0.25 


±0.009 


±0.016 


DMCyg 


223 


106 


54035.4065 


0.41986367 


0.523 


0.0066 B 


2.49 


0.091 


0.217 






(1900-2006) 


± 0.0009 


± 0.00000007 


+ 0.010 




±0.05 


±0.002 


±0.004 


V684 Cyg 


36 


53 


44253.2522 


0.54099474 


5.222 


0.0227 


19.30 


0.705 


1.304 






(1927-1980) 


± 0.0084 


±0.00000119 


+ 0.314 




± 1.16 


±0.042 


±0.079 


SUDra 


189 


103 


54111.5136 


0.66042294 


0.561 


0.0124 


1.70 


0.062 


0.094 






(1904-2007) 


± 0.0020 


± 0.00000019 


± 0.034 




± 0.10 


±0.004 


±0.006 


SVEri 


54 


104 


53998.8128 


0.71387846 


20.72 


0.0429 


58.06 


2.121 


2.972 






(1904-2006) 


± 0.0153 


± 0.00000142 


± 0.304 




±0.85 


±0.031 


±0.044 


VZHer 


203 


107 


53967.4819 


0.44033227 


0.678 


0.0207 


3.08 


0.112 


0.255 






(1899-2006) 


+ 0.0031 


± 0.00000018 


± 0.023 




± 0.10 


±0.004 


±0.009 


RR Leo 


298 


109 


54124.4177 


0.45240129 


1.869 


0.0072 


8.26 


0.302 


0.668 






(1898-2007) 


± 0.0009 


± 0.00000006 


± 0.009 




±0.04 


±0.001 


±0.003 


TV Leo 


27 


89 


52373.1085 


0.67286445 


3.271 


0.0307 


9.72 


0.355 


0.528 






(1913-2002) 


± 0.0184 


± 0.00000147 


± 0.284 




±0.84 


±0.031 


±0.046 


ULep 


33 


118 


54114.6303 


0.58147954 


0.685 


0.0077 


2.36 


0.086 


0.148 






(1899-2007) 


± 0.0019 


± 0.00000018 


± 0.030 




±0.10 


±0.004 


±0.006 


AVPeg 


443 


103 


54060.3926 


0.39038092 


0.894 


0.0084 


4.58 


0.167 


0.429 






(1903-2006) 


± 0.0010 


± 0.00000006 


± 0.010 




±0.05 


±0.002 


±0.005 


AR Per 


153 


107 


54124.4299 


0.42555066 


0.136 


0.0057 


0.64 


0.023 


0.055 






(1900-2007) 


± 0.0008 


± 0.00000006 


± 0.009 




±0.04 


±0.001 


±0.003 


RYPsc 


42 


96 


54037.5428 


0.52973431 


0.863 


0.0205 


3.26 


0.119 


0.225 






(1912-2006) 


± 0.0086 


± 0.00000062 


± 0.098 




±0.37 


±0.014 


±0.026 


CS Ser 


27 


102 


53530.4879 


0.52679833 


0.747 


0.0142 


2.84 


0.104 


0.197 






(1903-2005) 


± 0.0090 


± 0.00000050 


± 0.066 




±0.25 


±0.009 


±0.017 


BB Vir 


25 


94 


53849.6121 


0.47110683 


1.083 


0.0193 


4.60 


0.168 


0.357 






(1903-2005) 


± 0.0062 


± 0.00000060 


± 0.089 




±0.38 


±0.014 


±0.029 
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Table 5. Parabolic elements for RRab stars showing a well-defined linearly decreasing period. The note beside the value of the 
standard deviation indicates the stars having a known Blazhko effect and discussed in the text. 



Star 


N max 


Time coverage 


Epoch 


Period 


Quad, term 


s.d. 


dP/dt 




ct 






[years] 


[HJD-2400000] 


[d] 


[10" 10 -d/d] 


[d] 


[10" 10 -d/d] 


[dMyr 1 ] 


[Myr 1 ] 


SW And 


393 


111 


54093.3336 


0.44226187 


- 1.013 


0.0132 s 


-4.58 


-0.167 


-0.378 






(1894-2005) 


± 0.0016 


+ 0.00000010 


+ 0.013 




±0.06 


±0.002 


±0.005 


SX Aqr 


97 


93 


54018.3651 


0.53570932 


-0.537 


0.0139 


-2.00 


- 0.073 


-0.137 






(1913-2006) 


+ 0.0031 


± 0.00000025 


± 0.040 




± 0.15 


±0.005 


±0.010 


BR Aqr 


66 


105 


54048.3622 


0.48187046 


- 0.950 


0.0128 


-3.94 


-0.144 


- 0.299 






(1901-2006) 


± 0.0045 


± 0.00000024 


±0.031 




±0.13 


±0.005 


±0.010 


CP Aqr 


82 


72 


54024.3267 


0.46340227 


-0.516 


0.0103 


-2.23 


-0.081 


-0.176 






(1912-2006) 


± 0.0025 


+ 0.00000017 


+ 0.030 




±0.13 


±0.005 


±0.010 


TWBoo 


92 


92 


53918.4570 


0.53226977 


-0.503 


0.0095 


- 1.89 


-0.069 


-0.130 






(1914-2006) 


± 0.0015 


+ 0.00000019 


± 0.036 




±0.14 


±0.005 


±0.009 


AH Cam 


109 


56 


54119.3136 


0.36871531 


- 2.397 


0.0167" 


- 13.00 


- 0.475 


- 1.288 






(1951-2007) 


± 0.0029 


+ 0.00000039 


± 0.076 




±0.41 


±0.015 


±0.041 


APCnc 


24 


72 


52635.7481 


0.53293273 


-2.178 


0.0259 


-8.17 


- 0.299 


- 0.560 






(1930-2002) 


± 0.0192 


+ 0.00000138 


± 0.249 




±0.93 


±0.034 


±0.064 


WCVn 


110 


105 


54121.6211 


0.55175472 


-0.412 


0.0073 


- 1.49 


- 0.055 


- 0.099 






(1902-2006) 


± 0.0012 


± 0.00000014 


± 0.024 




±0.09 


±0.003 


±0.006 


RVCap 


258 


95 


51982.7919 


0.44773748 


- 1.260 


0.0270 s 


-5.63 


-0.206 


- 0.459 






(1906-2001) 


± 0.0144 


± 0.00000069 


± 0.084 




±0.37 


±0.014 


±0.030 


RXCet 


120 


111 


52172.1923 


0.57368560 


- 1.950 


0.0280 


-6.80 


- 0.248 


- 0.433 






(1890-2001) 


+ 0.0128 


± 0.00000061 


+ 0.075 




±0.26 


±0.010 


±0.017 


RZCet 


39 


77 


54036.5318 


0.51059118 


-4.032 


0.0336 


- 15.79 


-0.577 


- 1.130 






(1929-2006) 


+ 0.0133 


+ 0.00000127 


+ 0.224 




±0.88 


±0.032 


±0.063 


S Com 


89 


96 


54118.5568 


0.58658461 


- 0.773 


0.0097 


-2.64 


- 0.096 


-0.164 






(1911-2007) 


+ 0.0019 


+ 0.00000018 


± 0.032 




±0.11 


±0.004 


±0.007 


VXHer 


164 


90 


53919.4514 


0.45536088 


- 0.922 


0.0107 


-4.05 


-0.148 


-0.325 






(1916-2006) 


+ 0.0016 


± 0.00000013 


± 0.022 




±0.09 


±0.003 


±0.008 


V394 Her 


105 


77 


53932.4381 


0.43605116 


- 0.748 


0.0205 


-3.43 


-0.125 


- 0.287 






(1929-2001) 


+ 0.0049 


± 0.00000042 


+ 0.070 




±0.32 


±0.012 


±0.027 


RZLyr 


396 


111 


53983.3800 


0.51123785 


- 1.312 


0.0172 s 


-5.13 


-0.188 


-0.367 






(1895-2006) 


+ 0.0024 


± 0.00000014 


+ 0.019 




±0.08 


±0.005 


±0.003 


VI 095 Oph 


32 


64 


48832.4223 


0.45877735 


- 1.246 


0.0212 


-5.43 


-0.198 


- 0.433 






(1928-1992) 


± 0.0074 


± 0.00000091 


±0.174 




±0.76 


±0.028 


±0.060 


V964 Ori 


24 


104 


54080.5994 


0.50463562 


- 2.730 


0.0087 


- 10.82 


-0.395 


-0.783 






(1902-2006) 


+ 0.0045 


± 0.00000028 


± 0.039 




±0.15 


±0.006 


±0.011 


BHPeg 


131 


73 


54059.2925 


0.64098176 


-2.785 


0.0168 s 


-8.69 


-0.317 


- 0.495 






(1932-2006) 


+ 0.0028 


± 0.00000035 


± 0.084 




±0.26 


±0.010 


±0.015 


SW Psc 


30 


54 


39410.3459 


0.52124757 


-5.395 


0.0272 s 


- 20.70 


-0.756 


- 1.451 






(1912-1966) 


+ 0.0079 


± 0.00000152 


± 0.433 




± 1.66 


±0.061 


±0.116 


HKPup 


38 


76 


54100.6304 


0.73421141 


-7.216 


0.0285 


- 19.66 


-0.718 


- 0.978 






(1930-2006) 


+ 0.0161 


± 0.00000165 


±0.412 




± 1.12 


±0.041 


±0.056 


AT Vir 


25 


89 


53487.4580 


0.52577511 


-3.350 


0.0140 


- 12.74 


-0.466 


-0.885 






(1915-2004) 


+ 0.0049 


± 0.00000047 


± 0.086 




±0.33 


±0.012 


±0.023 



Sometimes, in the presence of large O-C variations we could 
not reach a uniquevocal evaluation of the elapsed cycles, and 
therefore we did not consider the star for the subsequent analy- 
sis. During this process we also corrected some times of maxi- 
mum brightness for obvious typographical errors; outliers (per- 
haps originated by errors in the decimal numbers) were also re- 
jected, especially within otherwise well-defined small groups of 
maxima. 

In several cases we corrected the linear elements since the re- 
spective O-C values were linearly increasing or decreasing. The 
new periods, calculated on a very long time baseline, provided 
more reliable O-C values. After this refinement, we started the 



analysis of the long-term behaviors. This procedure seems more 
reliable than comparing mean periods derived from groups of 
observations separated by very large gaps over a very long time 
baseline, as happens for globular cluster stars. Of course, in the 
case of galactic RR Lyr stars we met large gaps in the data, but 
we disregarded unclear cases. 

Table [2] summarises the inventory of the different patterns 
we observed in the O-C plots. Amongst the 123 cases, 81 RRab 
stars are brighter than 12.5 at minimum. To evaluate the com- 
pleteness of our sample, we note that the GCVS lists 164 RRab 
stars with a magnitude at minimum brighter than 12.5. Thus, 
our sample is complete at the 49% level; of course, there is no 
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way to recover the remaining 51%, since for these stars an in- 
sufficient number of times of maximum is available. Only in the 
future, with the progressive monitoring ensured by telescopes 
as TAROTs, will it be possible to increase the number of well- 
observed RRab stars. We note that amongst the 84 unstudied 
variables, 56 have 6 < 0°. 

The RRab stars showing a constant period for which we can 
provide improved elements are listed in Table[3] Formal errors at 
the 1 cr-level from the least-squares fits are reported. The RRab 
stars whose elements cannot fit all the observed times of max- 
ima are the most interesting cases. They are described in the next 
sections. We note that in spite of the short periods, RS Boo, AH 
Cam and EZ Cep are RRab stars. In particular, the GCVS clas- 
sification of EZ Cep (RRc) is wrong; the hipparcos photometry 
clearly shows a very steady ascending branch, typical of RRab 
stars (see Fig. 10 in Poretti |200T1 ). 

4. Stars with a linearly variable period 

The most striking result is the detection of a large number of 
O-C plots showing a parabolic pattern when linear elements 
are used. The cause of this pattern is a regular period variation. 
The coefficients of the parabolic fits and their error bars have 
been calculated by using standard routines (Press et al. 1 19921 1. 
Sometimes they are not well defined; as a rule of acceptance the 
correlation coefficient r of the least-squares parabola has to be 
larger than 0.70. The times of maximum have been fitted with 
the elements 

C = T + Po-E + a 3 -E 2 (1) 

where Pq is the period (in days) at epoch Tq and E is the cycle 
counter. Then, the linear rate dP/dt (in [d/d]) is derived from 

dP/dt = 2 ■ a 3 /{P) (2) 

where (P) is the average period, i.e., the ratio between the 
elapsed time (the difference between the first and last calculated 
maxima) and the elapsed cycles. In the literature it is more com- 
mon to indicate the rate of period changes using the f3 parameter, 
expressed in [d Myr 1 ] (e.g., Lee [19911 Jurcsik et al. 120011) . The 
ft values can be calculated from the second-order fits by apply- 
ing 

j3 = 0.0732- 10 10 -a 3 /(P) (3) 

It indicates how much the period will change in 10 6 years. In 
turn, this change can also be expressed as a fraction of the pul- 
sation period (e.g., Smith & Sandage [l9811 l, i.e.: 

a=pl(P) (4) 

Table [4] and Table [5] list the parabolic elements for RRab stars 
showing an increasing or decreasing period, respectively. With 
the aim of providing elements that will be directly useful in the 
planning of new observations, in each dataset the E=0 value has 
been assigned to the most recent time of maximum. Blazhko 
variables are also identified. The effect of having introduced the 
quadratic term can be evaluated for each star by comparing the 
lower and upper panels in Fig. [2] and Fig. [3] 

Amongst the 27 RRab stars showing a linearly increasing 
period (Tableland Fig. [2J, some cases are at the limit of ac- 
ceptance (0.70 < r < 0.76). In particular, the reliability of 
the older maxima is decisive in inferring the parabolic trends for 
V341 Aql, while that of the more recent ones is decisive for EZ 
Cep. In all the other cases the parabolic trend is clearly defined; 



in particular, r >0.85 in 20 cases. The standard deviations of the 
parabolic fits are less than 0.025 d in all cases except three (TV 
CrB, SV Eri and TV Leo). In particular, SV Eri is the RRab star 
showing the largest O-C variation and the largest dP/dt value, 
one order of magnitude larger than the average. 

Amongst the 21 RRab stars showing a linearly decreasing 
period (Table [5] and Fig.O the reliability of the older maxima is 
decisive in inferring the parabolic trends for TW Boo and V394 
Her. In the RV Cap case the parabolic trend is supported by both 
old and recent maxima, while the Blazhko effect (see next sec- 
tion) strongly contributes to reducing the r value (0.69). In 17 
cases we have r >0.85. 

The standard deviations listed in Tableland Table [5] have a 
mean value of 0.016 d and a median value of 0.014 d. They can 
be considered as an estimate of the accuracy of the photographic 
and visual times of maximum brightness, since they constitute 
the largest fraction of the data. We can infer that the weighting 
procedure would not have been very sensitive to the more accu- 
rate (by a factor from 5 to 10) CCD and photoelectric data; these 
types of data are still rare and grouped in the last few decades. 
This could explain why we obtained similar results by applying 
the weighting procedure when possible (see Sect. [3}. The above 
statistics and comparisons are quite satisfactory, and though the 
classification of some cases is somewhat doubtful, in general the 
observational scenario of stars with a linearly variable period is 
very well defined. 

Just a few of the stars listed in Table [4] and Table [5] were 
known to show a second-order term (BN Aqr, BR Aqr, CP Aqr, 
RW Cnc, W CVn, VZ Her, RR Leo, AV Peg, and AR Per), and 
in all cases there is an excellent agreement between the previous 
values (see for instance Olah & Szeidl 119781 and Szeidl et al. 
|1986l or, more in general, the GCVS, Kholopov [1988t and ours. 



4.1. Other types of O-C variations 

In our sample we can find examples of O-C patterns that are 
too complicated to be described with a polynomial fit. The O- 
C values derived from linear elements (not listed because they 
are valid only as rough average values on the considered time 
baseline) for several stars show irregular oscillations, with dif- 
ferent behaviors (Fig. 2}. We note that in several cases there is 
no gap in the observations and therefore there is no ambiguity 
about the general O-C behavior. The non-repeating O-C oscil- 
lation shapes rule out light-time effects. A Blazhko effect can 
be partially responsible for some scatter, but not for the general 
trend. 

The cases of SZ Hya, VV Peg, V759 Cyg and RU Cet are the 
simplest ones, since it seems that a single, abrupt change in the 
period has occurred. This also seems to be the case of UZ CVn; 
we rejected the parabolic fit since the left branch seems steeper 
than the right one, despite of the large scatter. XZ Cyg, RW Dra, 
RR Gem and RU CVn are a little more complicated; the changes 
occurred more than once, but still they are very sharp. The com- 
plexity of the O-C patterns increases in the cases of UY Boo, 
AR Her, Z CVn, ST Vir, AQ Lyr, RR Lyr itself, and RY Com, 
with several changes in the direction of the period variations. 
Other cases (EZ Lyr, RV UMa, BK Dra, SS Tau, BN Vul) show 
a smaller amplitude of the O-C values, resulting in more scat- 
tered plots. However, the general behaviors are very similar to 
the previous ones, with two or more small changes in the peri- 
ods. 
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5. Searching for Blazhko and light-time effects 

The most obvious way to verify the reliability of the new ele- 
ments, both linear and quadratic, is to analyze the residual O- 
C values. To do that, we applied the iterative sine-wave least- 
squares method (Vanicek l 19791 ). The frequency analysis was per- 
formed in the interval 0.00-0.20 cd~' . 

In several cases (SX Aqr, X Ari, RS Boo, RW Cnc, S Com, 
SU Dra, RZ Lyr, AV Peg, BH Peg and RY Psc) a very long pe- 
riod has been detected. Such a periodicity is often about half 
the time coverage and the amplitude is quite small. Therefore, 
we can infer that the parabolic fit is not perfect and that some 
systematic residuals were left, producing a spurious peak in the 
subsequent power spectrum. This feature can be noted looking 
at the O-C plots of the stars mentioned above in Figs. [2] and 
[3] the two branches seem to be asymmetric. We improved the 
O-C fitting by calculating a third-order least-squares fit, but 
we obtained a significant third-order term in three cases only 
(RS Boo, RW Cnc, SU Dra). They produce similar rates of the 
period changes: 1.70T0 -10 d/d (second-order fit) against 1.34 
•lO -10 d/d (third-order fit) for SU Dra, 7.08-10- 10 d/d (second- 
order fit) against 6.19 10" 10 d/d (third-order fit) for RW Cnc, 
1.42-10" 10 d/d (second-order fit) against 1.70 TO" 10 d/d (third- 
order fit) for RS Boo. We note that new TAROT observations in 
the next decade will allow us to disentangle the O-C behavior in 
these intriguing cases. 

RZ Cet seems to be a particular case. The residual O-C 
plot is quite unusual; a long term oscillation is clearly vis- 
ible and seems to be superimposed on the parabolic trend 
(Fig. [3}. The frequency analysis detects a low-frequency term 
at /=0.00008 cd _1 , corresponding to P= 12500 d, i.e., half the 
time coverage. The possibility of an erratic variation of the pe- 
riod is still plausible, but here the full amplitude is very large 
(0.082 d, i.e., 2.0 h). The standard deviation after the introduc- 
tion of the long period decreases from 0.034 d to 0.017 d and 
the fitting curve is perfectly sine shaped. These clues support the 
possibility of a light-time effect, but the gap in the O-C coverage 
hampers complete certainty. 

The presence of a Blazhko effect has been reported for SW 
Boo, SW Psc and FK Vul, but available times of maximum light 
are too scant to perform a reliable frequency analysis. The fre- 
quency analysis of the residual O-Cs related to SW And, DM 
Cyg, RZ Lyr and BH Peg did not detect any evident peak in 
the power spectra, though these RR Lyr variables are reported to 
be Blazhko stars. However, in some stars the amplitude of the 
Blazhko period changes over time, while in others the Blazhko 
effect might be more evident in the periodic modulation of the 
amplitude than in the O-C values. Therefore, in both cases the 
amplitudes of the O-C shifts could be too small to rise above 
the noise level. In other cases the detection has been positive. 
They constitute a convincing demonstration of the quality of the 
database. 

RV Cap - Kholopov (1988]) reports a Blazhko period 
of 225.5 d. Our frequency analysis refined this value to 
222.96+0.14 d. Figure [5] shows the power spectrum (lower 
panel) with the peak at f-0. 00448 cd~' clearly standing out. 
The full amplitude of the O-C variation is 0.045 d. In spite of the 
heterogeneous observational techniques, the O-C curve appears 
quite evident when folding the O-C over the 222.96 d period 
(upper panel; To=2430027.87). The Blazhko effect is responsi- 
ble for the large scatter observed in the O-C plot of RV Cap in 
Fig.[3]and for the low r value (0.69). 

AH Cam - The Blazhko period of AH Cam is one of the 
shortest known so far. Smith et al. d 19941 ) identified a period of 
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Fig. 5. The long-period Blazhko effect of RV Cap. Lower panel: power 
spectrum of the O-C values with respect to the parabolic elements; the 
highest peak is at 0.00448 cd~'. Upper panel: the same O-C values 
folded with the Blazhko period of 222.96 d. O-C measure unit is days. 



11 d by means of two photometric runs in 1989-90 and 1991— 
92 (covering 17 and 13 nights, respectively). The light curves 
show evident changes in shape and amplitude. The Blazhko ef- 
fect is also recoverable in the times of maximum light spanning 
56 years (Fig. [6}. The frequency analysis unambiguously identi- 
fied the peak at 0.0894 cd~' (lower panel), and the O-C plot over 
the refined period at 11.1808+0.0004 d (T =HJD 2438040.96) 
shows a full amplitude of 0.024 d (upper panel). 

RS Boo - The power spectrum reveals the low-frequency 
peak accounting for the slightly asymmetric O-C behavior (see 
above). However, a second peak is visible at 534.24 d, very close 
to the value given by Oosterhoff (119491 ), i.e., 533 d. The full 
amplitude is 0.010 d. On the other hand, there is no evidence of 
the additional period in the 58-62 d range suggested by Kanyo 
(119801 ); its amplitude, if real, is probably too small to be detected. 

TT Cnc - The power spectrum shows a peak at 
/=0.01145 cd _1 , corresponding to 87.35 d, close to the value 
determined by Szeidl ( 1 19681 ), i.e., 89 d. The resulting O-C curve 
shows some scatter. 

CX Lyr - The standard deviation of the O-C values with re- 
spect to the parabolic elements is quite high (0.036 d). The spec- 
tral window of the data is poor and hence the frequency analysis 
detects several peaks with comparable heights. The most con- 
vincing peaks correspond to periods of 227 d and 128 d. The 
O-C values folded over the two periods yield very similar plots, 
with a small preference for the 128 d curve; the full amplitude 
is 0.070 d and the standard deviation decreases to 0.027 d. No 
Blazhko effect was previously reported on CX Lyr. 

6. Period variations as evolutionary effects 

The dP/dt values listed in Tableland Table|5]constitute the first 
measurements of the rate of period changes in the sample of the 
best observed RR Lyr stars belonging to the Milky Way field. It 
is noteworthy that all the values, both positive and negative, have 
the same order of magnitude and therefore they should have a 
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Fig. 8. Examples of tracks showing blueward evolution at the begin- 
ning of the horizontal branch evolution. Triangles (in black) indicate 
a track with Z=0.0001 and M=0.86 M©, filled circles (in red) a track 
with Z=0.0004 and M=0.74 M©, empty circles (in blue) a track with 
Z=0.0007 and M=0.72 M©, squares (in green) a track with Z=0.0010 
and M=0.68 M©. The box indicates the instability strip domain (Lee 
[19931 
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Fig. 9. The distribution of the [Fe/H] values (as reported by Layden 
I1994t of the stars showing decreasing periods (16 out of 21 cases). 
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Fig. 7. Distribution of the rates of period changes. Constant periods 
(i.e., /?=0.0) are not shown. Note the extreme case of SV Eri in the 
positive part. 



common origin. Such an origin cannot be ascribed to a particu- 
lar cause (as it could for period variations induced by tidal ef- 
fects or by particular instabilities), but should be provoked by 
a quite general phenomenon, involving a whole population of 
stars. Though we cannot rule out the possibility that an individ- 
ual star amongst those listed in Tableland Table[5]is in a partic- 
ular transition phase (see Sect. 16.3b , the immediate answer is that 



the period variations are caused by long-term stellar evolution; 
indeed, RR Lyr stars change their radii and in turn their funda- 
mental radial periods when moving along the horizontal branch. 
The order of magnitude of the period changes we detected is that 
expected from evolutionary models: Lee dl99U reports ft values 
in the range 0.0-0.3 d Myr~' . Moreover, our values are similar to 
those reported by Jurcsik et al. ( 120011 ) in the case of a> Cen RR 
Lyr stars (J3 mean rate of 0.15 d Myr~') and discussed as evo- 
lutionary changes. Therefore, our analysis, providing a reliable 
estimate of the rate of period changes in field RR Lyr stars, al- 
lowed us to perfom a closer investigation of the stellar evolution 
along the horizontal branch. 
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6. 1 . Comparison between stars with period changes of 
opposite sign 

The histograms shown in Fig. Q summarise the properties of 
our sample. They provide evidence that a large fraction of RR 
Lyr stars experience appreciable changes in their structure in the 
time now covered by observations. The first immediate result is 
that the blueward evolution (decreasing periods) is as common 
as the redward one (increasing periods). Consequently, RR Lyr 
stars seem to be an excellent laboratory to verify the modifi- 
cations of the internal structure caused by stellar evolution un- 
der different conditions (e.g., contraction and expansion). On the 
other hand, we remind the reader that in 54 cases (i.e., a number 
of cases similar to the number of RR Lyr stars showing period 
variations) we did not detect a significant period variation. No 
period variation could occur at the turning point in the evolu- 
tionary tracks, but a constant period probably means that there 
are phases in the instability strip crossing when changes are not 
large enough to be measurable in 100 years or so. 

When limiting to cases where a period variability is noted, 
we get a similar distribution between negative and positive rates 
(Fig. [7} for our field RRab stars. We also note that the period 
ranges are nearly coincident in the two samples; in particular, 
we have a mean and median value of 0.53 d for stars with in- 
creasing period, and a mean and median value of 0.5 1 d for stars 
with decreasing period. This similarity is not observed for RRab 
stars in galactic clusters (Lee 119911 Jurcsik et al. |2001| l, where 
positive rates are largely more frequent. 

In our opinion, the large number of galactic stars showing 
decreasing periods is a strong confirmation of the evolution- 
ary tracks provided by Sweigart ( 1983 ar, d Lee & Demarque 
( 11990b . Most of these tracks cross the instability strip at high lu- 
minosities in the redward direction. However, Fig. [8] shows that 
for some combinations of masses and metallicities (e.g., by in- 
creasing the mass of the model with a lower metallicity), a blue- 
ward evolution at the beginning of the horizontal branch phase is 
predicted, both at high luminosity (Z=0.0001 and M=0.86 M ) 
and at low luminosity (Z=0.0010 and M=0.68 M Q ), and such 
blueward paths actually cross the instability strip domain. The 
theoretical sequences calculated by Dorman (1992) predict a 
blueward evolution at slightly lower luminosities, crossing the 
instability strip in fewer cases (e.g., Fig. 7 in Jurcsik et al. 1200 11 1. 
In general, all theoretical models predict a blueward evolution in 
connection with the prevalence of the H-shell burning over the 
He-central burning. The fact that the intersection between these 
blueward tracks and the instability strip is predicted only for par- 
ticular combinations of physical parameters is compensated by 
the longer lifetimes on the blueward path relative to the redward 
one (see next section). Therefore, once we have established that 
the blueward crossing occurs, the slow migration could explain 
why we observe RRab stars in this phase. 

To investigate the real agreement between the input param- 
eters of the evolutionary tracks shown in Fig. [8] and the stars 
of our sample, we considered the [Fe/H] values determined by 
Lay den (1994) by means of homogeneous low-to-moderate dis- 
persion spectra. Sixteen stars with decreasing periods are in- 
cluded in Layden's sample; we get [Fe/H]< -1.20 in 13 cases 
and a median value of [Fe/H]= -1.48 (Fig. [9]). These [Fe/H] val- 
ues t ransla te into Z < 0.0010 and Z = 0.0007 (e.g., Pietrinferni 
et al. 12004b . Therefore, the observed metallicities are in the same 
range of those used to generate the tracks shown in Fig. [8] A sim- 
ilar distribution is obtained for stars with increasing periods; in 
20 cases out of 24 we observe [Fe/H]< -1.20 and the median 
value is [Fe/H]= -1.52. This supports the hypothesis that all the 



stars with increasing and decreasing periods have similar char- 
acteristics, but the latter ones are less evolved than the former. 

6.2. The rate of period changes 

To avoid influence from a few very large positive or negative 
rates, we calculated the median values of the rates of the period 
changes to characterise the two distributions shown in Fig. [7] We 
obtained f3 = +0.14, a = +0.26 for the 27 stars with increasing 
periods and fi = -0.20, a = -0.43 for the 21 stars with decreas- 
ing periods. As already shown in early investigations by Sandage 
(11957b and then confirmed by Lee ( 11991b and Pietrinferni et al. 
(2004), the estimate of the horizontal branch lifetime is about 
~10 8 years, and a star spends -1/3 of its life in the instability 
strip. Moving uniformly through the strip with time, the corre- 
sponding rate of period change should be around a=0.02: Fig. [7] 
provides evidence of rates one order of magnitude greater. 

However, we have to take into account that when observing 
a star for 100 years we are monitoring a fraction correspond- 
ing to 10~ 5 of the lifetime spent crossing the strip. The rates 
of the period changes we measure are a sort of instantaneous 
value, which refers to a particular stage in the horizontal branch 
phase. Indeed, it is quite obvious that the stars cannot maintain a 
constant rate of some 0.1 d Myr -1 for several 10 6 years, since 
their periods would assume values largely outside the RRab 
range (0.4-0.8 d) and thus physically meaningless. Therefore, 
we should definitely refrain from considering the observed rates 
as the constant rate of the instability strip crossing of a specific 
stellar model. Such large variations (units of 0.1 d Myr~') must 
last for a few 10 6 years only. On the other hand, the 54 stars 
showing a constant period could be in the middle of their long 
crossing of the instability strip. 

The obvious interpretation of Fig|7] is to consider that the 
observed rates are those corresponding to the phases when the 
stars change their periods with the fastest rate. The rate spread is 
caused both by the different physical parameters of the stars and 
by the small differences in age. When considering the evolution- 
ary rates provided by Lee (11991b , it seems that large /? rates oc- 
cur only at the end of the instability strip crossing (i.e., roughly, 
1/20 of its horizontal branch lifetimes). However, such a propor- 
tion would imply a huge number of RR Lyr stars with nearly 
constant periods and a few stars with large (3 values. In our sam- 
ple we found a 2: 1 ratio. To explain the low number of stars with 
constant period we should admit that amongst the galactic RR 
Lyr stars there is a predominance of low-mass stars: the theoret- 
ical models predict (see Fig. 1 in Lee 1 199 il l that these stars cross 
the instability strip very quickly, close to the end of their hori- 
zontal branch lifetime. The case of S V Eri can support this point. 
It has a long period (0.71 d), close to the observational limit of 
RR Lyr periods. The theory predicts a large positive rate in this 
case and indeed we observe c=2.97, actually a huge value (see 
below). XX And has a similar period value (0.72 d), but its rate 
of period change is one order of magnitude smaller (c=0.20). 
Thus, SV Eri seems to have entered the rapid final phase, and it 
should be considered a very peculiar object experiencing a rapid 
final evolution. 

The median values also suggest that galactic RR Lyr stars 
experienced slightly faster period decreases than increases. To 
evaluate the physical meaning of this fact, we need to com- 
pare them more closely with the lifetimes of the blueward and 
redward crossings. Actually, the blueward crossings are about 
1 .6 times longer than the redward. As an example, the two tracks 
with Z=0.0004 and Z=0.0007 in Fig. [8] (i.e., those indicated by 
filled and empty circles, respectively) yield 33 and 28 Myr for 
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the blueward evolution and 21 and 17 Myr for the redward evo- 
lution. Moreover, the stars should change their radius to a lesser 
extent in the blueward evolution than in the redward one (0.13 
against 0.20, in log R). Therefore, at first glance, the positive /3 
values should be on the average greater (also in absolute value) 
than the negative ones, since a larger period variation occurs in 
a shorter time. This expectation seems to be contradicted by the 
histograms shown in Fig. [7] The contradiction could be partially 
explained by taking into account the observation that the short 
lifetimes on the final part of the redward evolution strongly re- 
duce the possibility of observing these stars. Nevertheless, the 
RRab stars showing rapid period decreases still need more de- 
tailed theoretical investigations. For example, their large number 
could result from a burst in the stellar formation process in our 
galaxy. 

Finally, we have to address a possible discrepancy between 
the observed and predicted rates of the period changes. The the- 
oretical evolutionary rates calculated by Lee d 1 99 lb and recalled 
above predict /3 values smaller than 0.30 for almost all the stars. /3 
values larger than 0.30 (stars with very low mass, M = 0.68 Mq, 
and well-defined metallicity, Z=0.0004) are only reported in ex- 
ceptional cases and the lifetimes of these phases are very short. 
In our sample we have 18 stars with < f3 < 0.30 and 9 stars 
with f3 > 0.30. Also in this case, the 2:1 ratio seems too small 
to be caused by the presence of exceptional cases. Furthermore, 
the a-2.91 value of SV Eri recalled above is much higher than 
expected. Smith & Sandage (119811 also quoting Iben 1 197TT > in- 
dicate a-Q.5 as an extreme value, but such a value is not an 
uncommon one (see Hig. [7j an d is much smaller than the SV Eri 
value. 

Our conclusion is that the theoretical models should be fine 
tuned to match the observed rates in a better way, now that the 
extensive GEOS database allows us to measure them in a reliable 
sample. 

6.3. Erratic changes 

During the final phase of the helium core burning it is possible 
that instabilities occur and that period variations are not regu- 
lar, both around a constant value or superimposed on the gen- 
eral trend (Sweigart & Demarque 119731 ). This possibility has 
been suggested by Sweigart & Renzini (1 19791 1 and has been in- 
voked by Lee ( 119911 ) to explain the few cases of negative values 
in the rate of period changes. RR Lyr stars showing blueward 
evolution have been considered as suspicious, invoking period 
changes caused by temporary instabilities, but we are providing 
evidence that the number of RR Lyr stars subject to blueward 
evolution is similar to the number of stars subject to the red- 
ward evolution. Therefore, the erratic changes should be consid- 
ered as representative of a particular evolutionary stage, as the 
final helium core burning could be, but not as an ad hoc solu- 
tion to explain negative rates only. We also stress that the erratic 
changes occur in both directions (see Fig.|4|. In our opinion, the 
sudden, sharp period changes occurring in the cases similar to 
SZ Hya could be the result of particular temporary instabilities, 
which can also produce sudden blueward excursions (Sweigart 
& Renzini fl 9791 1. We also stress that only a comprehensive anal- 
ysis performed on a large database can avoid overlooking these 
sudden changes as evolutionary ones. 

7. Conclusions 

The analysis of O-C variations over a timescale of 100 years or 
more has proved to be a powerful tool for providing quantitative 



tests of the stellar evolution of the horizontal branch stars. We 
can stress some well-established observational facts: 

1. RRab stars experiencing blueward evolution (i.e., period de- 
creases) are quite common, only slightly less so than RRab 
stars experiencing redward evolution (i.e., period increases). 
The period ranges covered by the two groups are very simi- 
lar and the mean and median period values are nearly coin- 
cident; 

2. the absolute values of period changes are larger than ex- 
pected. Also in the extreme case of the rapid-evolving star 
SV Eri the rate is much larger than expected; 

3. the O-C behavior can be very complicated in some cases, 
with abrupt or irregular period variations, rather than mono- 
tonic. The regular variations caused by the light-time effect 
(and hence duplicity), often invoked to explain large O-C 
excursions, are not convincingly observed in our extensive 
sample. The physical explanations should be searched in the 
stellar structure; 

4. Blazhko effect is often superimposed on secular changes, but 
the monotonic trend due to evolutionary variations still re- 
mains visible. 

As a general conclusion about the comparison between our ob- 
servational results and the theoretical predictions, we claim that 
there is a very powerful feedback between the two approaches. 
In particular, theoretical investigations should take into account 
that we have observational evidence of many RR Lyr stars 
showing blueward evolution. The theoretical models should also 
match the observed j3 values in a more satisfactory way, as these 
seem to be higher than expected, both for redward and blueward 
evolutions. On the other hand, the observational effort to monitor 
RR Lyr stars should be continued, possibly extended to stars at 
fainter magnitudes, and the accuracy of the maximum time deter- 
minations should be improved (by using automated telescopes), 
thus obtaining the same information on a shorter time baseline. 
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Fig. 2. O-C values for RRab stars. Filled circles indicate O-C values calculated from the parabolic fit (period increasing at a constant rate), crosses 
from the linear fit (constant period). On the x— axis we plot the elapsed cycles (E/10000); on the ordinate one we plot the O-C values in days. The 
cases of V341 Aql, RS Boo, RW Cnc, EZ Cep and SU Dra are discussed in the text. 
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Fig. 3. O-C values for RRab stars. Filled circles indicate O-C values calculated from the parabolic fit (period decreasing at a constant rate), 
crosses from the linear fit (constant period). On the x— axis we plot the elapsed cycles (E/10000); on the ordinates the O-C values are in days. 
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Fig. 4. O-C values for RRab stars calculated with a linear fit and showing irregular variations. They are ordered by increasing complexity of their 
O-C patterns. On the x-axis we plot the elapsed cycles (E/10000); on the ordinate one we plot the O-C values in days. 



